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SCALE  SIZES  AND  LIFETIMES  OF  F  REGION  PLASMA  CLOUD  STRIATIONS 
AS  DETERMINED  BY  THE  CONDITION  OF  MARGINAL  STABILITY 

1.  Introduction 

Phenomenologists  have  been  pressed  to  account  for  certain  features  of 
late  time  barium  cloud  striations,  despite  the  early  time  success  and  near 
universal  acceptance  of  results  from  the  E  x  B  gradient  drift  instability  theory 
[Linson  and  Workman.  1970;  Perkins  et  al.,  1973,  Francis  and  Perkins,  1975]. 

We  refer  to  the  simultaneous  occurrence  of  three  phenomena  as  documented 
in  the  1977  STRESS  program:  a)  the  generation  of  structure  on  scales  as 
small  as  15  meters  [Baker  and  Ulwick,  1978] ;b)  the  apparent  ubiquity  of 
visible  striations  separated  by  distances  approaching  one  kilometer;  and 
c)  the  survival  for  hours  (sometimes  referred  to  as  "freezing  up")  of  this 
kilometer-scale  structure  as  evidenced  by  propagation  studies  [Prettie  et  al. , 
1977].  If  the  gradient  drift  instability  is  at  work  producing  small  struc¬ 
ture  through  a  series  of  bifurcations,  why  do  the  kilometer  scale  structures 
persist?  What  parameters  select  out  the  kilometer  scale  size  so  often? 

(Observers  have  often  noted  that  visible  structuring  seems  to 
halt  when  scale  sizes  transverse  to  the  neutral  wind  direction  decrease  to 
just  under  one  kilometer  [J.  A.  Fedder,  W.  Chestnut,  private  communication, 

1980 ]  .  Past  this  point,  there  is  a  tendency  for  the  striations  to  drift 

in  unison  as  long  as  they  can  be  seen.) 

We  support  the  gradient  drift  instability  as  the  explanation  for  stri- 

ation  behavior.  In  fact,  the  emergence  at  late  times  of  density  power  law 
spectra  from  numerical  simulations  [Scannapieco,  et  al.  1976]  is  in  agreement 
with  observation  [Baker  and  Ulwick,  1978].  We  feel,  however,  that  attempts 
to  explain  the  above  mentioned  three  features  have  been  hindered  by  approach¬ 
ing  the  problem  from  two  opposite  extremes:  1)  application  of  analytic 
Manuscript  submitted  September  29,  1980. 
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results  from  an  idealized  one  dimensional  slab  model;  and  2)  brute  force  case- 
by-case  numerical  simulation  beginning  with  an  initially  smooth  plasma  dis¬ 
tribution,  and  proceeding  toward  a  highly  structured  state.  One  is  not 
surprised  that  some  disagreement  exists  between  an  idealized  linear  result 
and  the  highly  structured  nonlinear  environment.  On  the  other  extreme, 
global  numerical  simulation  on  present  day  computers  must  eventually  be 
hindered  by  inadequate  resolution.  We  feel,  too,  that  the  role  of  diffusion 
(both  classical  and  turbulent)  has  not  been  sufficiently  emphasized  in  de¬ 
termining  late  time  structure.  This  is  partly  because  present  day  computers 
cannot  resolve  both  the  global  cloud  scale  and  the  diffusion  scale  due  to 
storage  and  speed  limitations.  In  order  to  fill  in  these  gaps,  we  shall  ex¬ 
tract  from  the  basic  one  level  (F  region)  striation  model  a  scaling  law 
which  allows  individual  structures  to  provide  information  concerning  "bifur¬ 
cation  tendencies"  applicable  to  other  similar  structures  of  arbitrary  size. 

We  demonstrate  in  Section  2  that  bifurcation  or  non-bifurcation  of  a  par¬ 
ticular  structure  hinges  upon  whether  a  diffusion  parameter  R,  analogous 
to  the  Reynolds  number  for  neutral  flows,  exceeds  a  critical  value.  This 
R  value  provides  an  estimate  for  the  scale  size  of  the  state  of  marginal  bi¬ 
furcation  tendency.  It  also  provides  an  estimate  for  the  lifetime  of  the 
marginal  state.  We  determine  the  critical  R  by  high  resolution  numerical 
simulation  of  an  isolated  structure  of  relevant  geometry.  Just  enough  physical 
diffusion  is  added  to  prevent  bifurcation  and  loss  of  resolution. 

We  are  aware  that  mechanisms  other  than  diffusion  can  affect  structure 
on  suitable  scales.  Among  these  are  inertial  forces,  kinetic  (finite 
gyroradius)  effects,  and  E-to-F  layer  coupling.  It  is  our  desire,  however, 
to  present  the  roost  basic  model  which  seems  to  agree  with  the  stated  obser¬ 
vations.  In  order  to  explain  the  side-by-side  existence  at  late  times  of 
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15  meter  and  kilometer  scale  sizes  [Baker  and  Ulwick,  1978]  we  hypothesize 
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that  the  cloud  Is  dominated  by  turbulent  diffusivity  of  order  100  m  /sec, 
but  that  selected  regions  may  reflect  non  turbulent  electron  diffusion  of 
order  1  m^/sec. 

In  Section  3  we  describe  the  numerical  techniques  by  which  the  equations 
of  motion  are  solved.  These  techniques  reflect  appropriate  advances  in  the 
state-of-the-art  since  the  early  simulations  of  large  scale  striation  morphology 
Lzabusky,  et  al. ,  19-73;  Scannapieco,  et  al.,  1974,  1976].  The  scaling 

derived  in  Section  2  allows  us  to  simulate  a  small  region  of  the  plasma  with 

adequate  resolution,  providing  that  boundary  conditions  can  be  specified 

appropriately.  Section  4  contains  results  of  the  numerical  simulations  and 

illustrates  how  the  critical  R  depends  upon  the  conductivity  ratio  M.  In 

Section  5  we  show  that  our  calculated  R  values  are  compatible  with  rocket 

probe  data  on  scales  of  tens  of  meters  [.Baker  and  Ulwick,  1978],  We  show 

that  if  turbulent  diffusivity  is  operative,  a  wide  range  of  conditions  will 

lead  to  marginally  stable  scale  sizes  of  approximately  one  kilometer.  In 

addition,  we  predict  that  the  lifetime  of  the  kilometer-scale  structure 
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should  be  of  order  10  sec.,  in  agreement  with  propagation  experiments 
conducted  during  the  STRESS  program  LPrettie  et  al.,  1977].  Thus  it  is 
possible  that  our  results  may  help  to  explain  the  "freezing  up"  of  visible 
kilometer-scale  barium  cloud  structures  [J.  A.  Fedder,  W.  Chestnut,  private 
communication,  1980]. 

2.  Equations  of  Motion  for  the  One  Level  Model 
The  two  dimensional  field  line  integrated  Pedersen  conductivity  model 
appropriate  to  F-region  clouds  consists  of  the  following  equations  (in  cgs 
units)  which  are  cast  in  a  frame  drifting  with  the  ambient  plasma: 

-■-7  .  EV  +  K7^£ 
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(1) 


V  <J>  x  z 


(2) 


V 


£ 

B 


V 


Z  7$  =  E 


(3) 


where  Z  is  the  magnetic  field-line -integrated  Pedersen  conductivity,  K  is 
the  cross-field  diffusivity  of  the  cloud  plasma,  V  is  the  local  plasma  drift 

relative  to  the  ambient  drift  velocity,  c  is  the  speed  of  light,  B  =  B  £  is 
the  constant  magnetic  field  strength,  4>  is  the  induced  electrostatic  potential, 
and  E  =  E  y  is  the  ambient  electric  field  in  the  rest  frame  of  the  neutral 
atmosphere.  Equations  (1)  -  (3)  are  two  dimensional  (x,y)  where  x  is  parallel 
to  the  relative  neutral  wind  and  y  is  the  direction  of  E^.  Except  for  the 
assumption  of  constant  K,  these  equations  are  accurate  to  first  order  in  the 
ratio  of  ion  collision  frequency  to  gyro frequency.  Equation  (1)  results  from 
multiplying  the  ion  continuity  equation  by  the  Pedersen  mobility  and  integrat¬ 
ing  along  the  magnetic  field.  One  assumes  that  the  electrostatic  potential 
and  therefore  V  are  constant  on  field  lines.  The  diffusion  term  in  (1)  results 
from  the  assumption  of  a  constant  Pedersen  mobility  and  K  value  within  the 
cloud.  The  assumption  of  constant  K  is  admittedly  simplistic.  For  a  typical 
barium  cloud  without  turbulence,  K  is  roughly  proportional  to  the  local  plasma 
density  [Perkins  et  al.,  1973].  For  a  turbulent  cloud,  K  may  be  scale  size 
dependent  [Goldman  and  Sperling,  1979].  Attempts  to  model  the  variation  of 
K  introduce  an  additional  scaling  parameter  and  tend  to  muddy  the  results 
without  changing  them  significantly.  Equation  (3)  is  the  equation  of  quasi¬ 
neutrality,  which  states  that  the  field  line  integrated  electric  current, 

J  =  Z  E  *  ~  7$)*  must  be  divergence  free. 

For  the  sake  of  completeness,  we  should  mention  that  a  two  layer  model 
with  one  layer  for  the  cloud  and  another  for  the  background  ionosphere  has 
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been  in  use  for  several  years  [Lloyd  and  Haerendel,  19733.  One  notices  that 
simulations  performed  with  the  model  tend  to  yield  late  time  configurations 
in  which  there  is  a  strong  correlation  between  density  distributions  in  the 
two  layers  [Scannapieco  et  al.,  1976;  S.  T.  Zalesak,  private  communication, 
1980].  This  may  lend  support  to  one  level  predictions  concerning  bifurcation 

of  a  particular  structure,  providing  that  the  conductivity  is  dominated  by 
the  Pedersen  component  at  the  cloud  level.  Equations  (1)  -  (3)  may  be  put 

into  dimensionless  form  as  follows  [McDonald,  et  al.,  1978].  Let 

x  =  L  xy 
—  o  — 

dt  =  t  dt' 
o 

e  -  eo  r 

I  -  V„  S'  (4) 

<f>  =  LE  <j>' 

o  o 

Vo  -  CIVBI' 

and  t  -  L  /V  • 

ana  o  o  o 

Here  L  is  a  measure  of  the  cloud's  gradient  scale  size,  and  all  primed 
o 

quantities  are  dimensionless.  V  is  the  relative  drift  speed  between  the 

o 

ambient  plasma  and  the  neutral  atmosphere,  and  is  the  ambient  integrated 
Pedersen  conductivity.  In  this  caper  we  take 


h  -  </*  dx  dy/'Vc  fez)2  dxdy)^ 


where  the  region  of  integration  is  a  rectangular  box  containing  the  structure 
of  interest.  We  have  experimented  with  more  sophisticated  definitions  of  Lq 
which  are  insensitive  to  the  size  of  the  box.  These  yield  results  within  20% 
of  those  given  by  (5)  for  the  cases  presented  here.  We  prefer  to  retain  the 


straightforward  simplicity  of  (5). 
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Substituting  (4)  into  (1)  —  (3)  and  dropping  primes  from  all  variables 
we  have  the  dimensionless  equations 


-^5  =  -  V.  I  V  +  R~^  £ 

(6) 

V  =  -  V  if)  x  z 

(7) 

V  .  EV  *  =  i  .  7  E, 

(8) 

where 

R  =  V  L  /K. 
o  o 

(9) 

Note  that  R  is  analogous  to  the  Reynolds  number  for  neutral  flows,  with  the 
kinematic  viscosity  replaced  by  the  cross-field  plasma  diffusivity. 

Equations  ( 6)  -  (8)  reveal  that  the  evolution  of  a  plasma  cloud  is 

completely  determined  by  initial  cloud  geometry,  boundary  conditions,  and  the 
value  of  R.  Thus  the  answer  to  whether  or  not  a  given  structure  will 
bifurcate  depends  upon  whether  R  is  greater  or  smaller  than  some  critical 
R  value  applicable  to  that  structure.  For  sufficiently  small  R,  diffusion 
will  be  dominant  and  will  smooth  out  small  structures  faster  than  they  can 
be  created.  For  sufficiently  large  R,  diffusion  will  be  negligible,  allowing 
steepening  and  bifurcation  to  proceed. 

If  the  critical  R  for  a  particular  structure  is  known,  we  can  use  (9) 
to  estimate  the  scale  size  of  the  marginally  stable  state,  providing  esti¬ 
mates  for  V  and  K  are  available: 


The  lifetime  of  the  marginal  state  may  also  be  estimated.  The  time  re¬ 
quired  for  diffusion  to  degrade  the  structure  is  approximately 


=  R  L  /V 
o  o 


(ID 


One  sees  that  tD  is  R  times  the  one  dimensional  gradient  drift  instability 
growth  time.  Anticipating  that  R  is  a  large  number,  (11)  implies  that  the 
marginal  state  can  persist  for  times  much  greater  than  basic  gradient  drift 
structuring  times. 

In  this  work  we  investigate  bifurcations  which  originate  near  the  tips 
of  striations.  This  seems  to  be  the  most  likely  place  for  new  structure 
to  emerge,  although  there  is  photographic  evidence,  for  example  in  SPRUCE, 
that  some  new  structure  can  emerge  from  sides  of  striations.  The  procedure 
for  estimating  the  critical  R  for  a  given  Initial  condition  on  Z,  will  be  to 
carry  out  a  set  of  simulations  from  equations  (1)  -  (3)  with  assorted  values 
of  K.  An  upper  limit  on  the  critical  R  will  be  taken  to  be  the  smallest  R 
at  the  time  of  f^rst  bifurcation.  A  lower  limit  will  be  taken  as  the  largest 
R  among  non-bifurcating  cases  at  times  used  to  determine  the  upper  limit. 

The  initial  condition  on  Z  will  be  taken  as  an  appropriate  form  representing 
a  single  .ated  striation  tip.  This  allows  numerical  resolution  to  be 

4  to  8  times  that  of  our  earlier  study  [McDonald  et  al.,  1978 J. 


3.  Numerical  Simulations 

Equations  (1)  -  (3)  are  advanced  in  time  using  numerical,  techniques 
described  below.  A  more  detailed  description  of  the  NRL  one-level  stria¬ 
tions  code  has  been  given  by  McDonald  et  al..  [1979],  We  choose  the  dimension¬ 
al  set  rather  than  (6)  -  (8)  to  facilitate  comparison  with  experiment. 
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With  £(x,y)  given  at  a  particular  time,  we  calculate  <J>  from  (3),  and  use  the 
result  in  (2)  to  arrive  at  the  flow  field  V.  Then  we  use  (1)  to  advance  £ 
by  one  timestep,  after  which  the  cycle  can  be  reinitiated.  Finite  difference 
representations  are  used  for  all  derivatives  in  (1)  -  (3). 

Solution  of  the  variable  coefficient  elliptic  equation  (3)  accounts  for 
nearly  two-thirds  of  simulation  time  for  grids  as  large  as  that  used  here 
(162  by  82  points).  We  use  the  Chebychev  semi-iterative  method  to  calculate 
a  series  of  approximations  which  converge  to  the  true  solution  of  (3) 

[Varga,  1962;  McDonald,  1980].  This  method  was  chosen  because  of  its 
efficient  execution  on  a  vector  computer  such  as  the  NRL  Texas  Instruments 
ASC.  It  also  allows  convenient  implementation  of  a  wide  variety  of  boundary 
conditions.  Derivatives  are  represented  bv  second  order  finite  differences 
on  a  uniform,  non-staggered  grid  (£,  <|>,  and  V  are  calculated  on  the  same 
set  of  points).  It  is  possible  to  obtain  an  exact  solution  to  the  finite 
difference  analog  of  (3)  [Madala,  1978],  but  computer  time  and  storage  re¬ 
quirements  become  excessive  for  large  grids.  We  use  time  extrapolation  to 

obtain  an  accurate  initial  approximation  to  <j>  ,  and  carry  out  iterations 

-4 

until  the  root  mean  square  residual  error  is  approximately  3  x  10  of  the 
root  mean  square  source  term.  Reduction  of  the  error  much  below  this  level 
would  not  necessarily  improve  results,  since  discretization  introduces  errors 
of  this  order  or  greater  into  the  equations. 

For  the  advection  term  in  equation  (j)  we  use  the  two  dimensional  flux 
correction  method  of  Zalesak  [1979].  This  technique  can  be  incorporated  into 
a  number  of  standard  advection  schemes  to  prevent  formation  of  spurious 
oscillations  in  the  advected  quantity.  There  are  two  reasons  why  these 
oscillations  should  be  avoided:  (1)  in  a  divergenceless  flow  field  (see 
eq.  (2))  new  maxima  or  minima  cannot  be  created:  and  (2)  high  order  schemes 
can  drive  £  negative,  resulting  in  degeneracy  of  eq.  (3)  and  physically 
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meaningless  velocities.  Traditional  remedies  include  local  "fill-in"  and 
addition  of  artificial  diffusion  globally.  The  flux  correction  technique 
operates  locally  in  regions  where  the  oscillations  tend  to  form,  effectively 
truncating  the  scheme  to  first  order  such  that  ripples  cannot  form.  The 
diffusion  term  in  equation  (1)  is  added  in  a  separate  step  using  second  order 
spatial,  forward  time  differencing.  In  the  integration  of  equation  (1)  we 
tacitly  subtract  from  V  the  structure's  centroid  velocity.  This  keeps  the 
interesting  structure  from  drifting  off  the  computational  grid.  More 
important,  it  reduces  the  number  of  grid  cells  through  which  the  structure 
drifts  during  a  simulation.  This  helps  to  reduce  discretization  error  due 
to  the  irreversible  loss  of  information  tfiat  occurs  when  an  arbitrary  pro¬ 
file  is  advected  at  constant  velocity  through  an  Eulerian  grid. 

The  simulations  were  carried  out  on  a  grid  of  162  by  82  points  in  the 
x  and  y  directions  respectively.  Each  simulation  was  carried  out  with  a 
uniform  grid  interval  of  80  meters  and  was  then  rerun  with  a  grid  interval 
of  40  meters  as  a  check  on  the  sensitivity  of  results  to  resolution.  The 
use  of  a  non-stretched  grid  increases  execution  efficiency  of  the  code.  The 
problems  arising  from  electrostatic  images  on  the  boundaries  are  partially 
alleviated  bv  the  use  of  transraittive  boundarv  conditions  described  below. 

The  calculation  procedure  Is  to  update  quantities  on  the  interior  of  the  grid 
(160  x  80  points)  by  the  methods  described  above,  and  then  calculate  ex¬ 
terior  values  from  the  boundary  conditions. 

Boundary  conditions  used  during  the  simulation  are  as  follows. 

Neumann  boundaries 

|-  =  0  (12) 
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involved  in  updating  equation  (l).  This  allows  inflow  and  outflow  of  ambient 
plasma  on  three  sides  of  the  grid  where  structure  does  not  intersect  the 
boundary  (see  Figures  1-4).  It  is  also  a  good  approximation  for  the  fourth 
side  (the  left  boundary  of  each  figure)  since  I!  contours  tend  to  align  with 
the  neutral  wind  direction  (x)  except  near  striation  tips.  For  the  present 
study,  the  Neumann  condition  seems  more  realistic  than  the  periodic  bounda¬ 
ries  used  in  most  earlier  simulations.  The  "re  entry"  of  plasma  resulting 
from  periodicity  would  be  nonphysical  in  the  study  of  an  isolated  structure. 

The  physically  proper  boundary  condition  on  4>  is  that  its  gradient  van¬ 
ish  at  infinity.  Since  we  are  constrained  to  a  grid  of  finite  extent,  we  can 
either  stretch  the  grid  to  move  the  boundary  far  from  the  structure  or 
attempt  to  match  the  interior  solution  with  an  appropriate  exterior  solution. 
We  have  chosen  an  approximation  to  the  second  alternative  and  have  developed 
a  physically  motivated  algorithm  which  is  admittedly  simplistic,  but  has 
been  validated  empirically.  Our  approximate  boundary  condition  is 


+  s  a_ 

dn  Sn 


=  0 


(13) 


where  n  is  the  outward  normal  and  S  is  a  constant  taken  to  be  the  cloud's 
initial  scale  size  (see  (15)).  This  states  that  the  normal  component  of 
the  polarization  electric  field  must  decrease  in  the  outward  direction  on 
a  spatial  scale  comparable  to  that  of  the  gross  structure.  We  choose  to 
cast  (13)  in  terms  of  the  first  two  derivatives  rather  than  $  and  its  first 
derivative  so  that  the  addition  of  a  constant  to  <£  would  have  no  effect 
upon  structure  evolution.  We  have  validated  (13)  by  comparing  analytic 
and  finite  difference  solutions  to  (3)  for  elliptic  "waterbag"  distributions 
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Z  (x,y).  The  condition  (13)  gave  results  superior  to  the  Utual  periodic, 
Neumann,  or  Dirichlet  boundaries. 

Initial  conditions  for  the  simulations  are  taken  to  be  descriptive 
of  the  tip  of  a  striation  (see  upper  portions  of  Figures  1-4).  Taking 
the  origin  at  the  center  of  the  rectangular  grid,  we  have  at  t  =  0 


1=1  +  (M-l)  exp(-y2/s2),  x  <  0 
=  1  +  (M-l)  exp(^(x2+y2)/S2  ,  x  ^  0 


(14) 


where  M  is  the  ratio  of  peak  conductivity  to  ambient.  For  all  cases 
presented  here, 

S  =  1  km  .  (15) 


The  computational  domain  for  results  shown  in  Figures  1-4  is  6.4  by  12.8 

kilometers.  That  we  have  normalized  Z  to  the  ambient  Z  has  no  effect 

o 

on  the  model  (1)  -  (3)  (i.e.,  Z  can  be  replaced  with  Z/Z  ).  All  other  varia- 

o 

bles  retain  their  dimensionality.  For  all  cases  we  take  B  =  .5  gauss,  and 
Eq  =  -1.67  x  10  2  stat  volts/cm,  so  that  Vq  =  100  m/s.  In  this  particular 
simulation  it  is  not  necessary  to  perturb  the  initial  condition,  since  (14) 
does  not  represent  a  steady  state.  One  must  note,  however,  that  discreti¬ 
zation  and  roundoff  errors  provide  effective  pertubations  to  any  simulation. 

4.  Results 

Simulation  results  used  to  determine  upper  and  lower  limits  on  R  for 
each  of  four  M  values  (2,  5,  10,  and  30)  are  shown  in  Figures  1-4.  For 
each  M,  a  set  of  simulations  were  carried  out  for  various  K  values  in  order 
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to  help  locate  the  demarcation  between  bifurcating  and  non  bifurcating  states. 
For  each  M  value  we  present  the  two  simulations  which  most  accurately  bound 
the  critical  R.  In  each  case  this  involves  a  pair  of  K  values  separated  by 
a  factor  of  2.  The  range  could  be  narrowed  by  further  simulations,  but  the 
approximations  involved  in  taking  K  constant  and  specifying  boundary  con¬ 
ditions  at  finite  distances  would  lend  some  uncertainty  to  the  result.  Our 
earlier  critical  R  value  for  an  M  =  11  structure  [McDonald  et  al. ,  1978]  is 
approximately  a  factor  of  2  lower  than  one  would  expect  from  the  present  work 
(see  Figure  6).  We  attribute  this  to  the  use  of  different  initial  conditions, 
boundary  conditions,  lower  spatial  resolution,  and  a  different  procedure 
for  determining  the  critical  R. 

In  each  of  Figures  1-4  the  left  column  gives  contour  plots  of  £ 
for  cases  where  K  is  small  enough  to  permit  bifurcation.  The  plot  just 
below  the  initial  condition  plot  shows  £  at  a  time  when  it  is  qualitatively 
clear  that  bifurcation  is  about  to  occur.  This  time  is  chosen  as  appro¬ 
priate  to  determine  the  critical  R  from  (5)  and  (9).  We  choose  not  to 
determine  R  at  t  =  0  because  the  uniform  separation  of  contours  is  not 
"natural."  This  becomes  clear  when  one  examines  the  non  bifurcating 
states  in  the  right  column  of  Figures  1-4.  A  quasi  stationary  state  is 
established  only  after  a  sufficient  amount  of  steepening  has  occurred.  It 
is  the  ftability  of  a  pre-steepened  and  thus  more  "natural"  state  that  we 
wish  to  address.  When  the  structure  is  near  bifurcation,  the  selection 
of  a  time  for  calculating  R  is  not  crucial.  As  shown  in  Figure  5,  the 
length  scale  tends  to  level  off  just  prior  to  the  development  of  new 
structure.  Thus  our  resalts  are  not  sensitive  to  the  precise  definition  of 
an  onset  time. 
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Fig.  1  —  Contour  plots  of  I(x,y)  at  selected  times  illustrating  the  demarcation  between  bifur¬ 
cating  (left  column)  and  non-bifurcating  (right  column)  states  for  M  ■  2.  Contours  are  spaced 
linearly  between  2  =  1  (ambient)  and  M.  Limits  on  the  critical  R  are  evaluated  at  times  when 
bifurcation  of  the  low  K  striation  (second  plot  in  left  column  for  Figures  1—4)  is  eminent. 
Times  appear  in  the  upper  right  of  each  plot. 
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Fig.  2  —  Contour  plots  of  I(x,y)  at  selected  times  illustrating  the  demarcation  between  bifur¬ 
cating  (left  column)  and  non-bifurcating  (right  column)  states  for  M  ■  5.  The  diffusion  dominated 
case  is  plotted  at  300  seconds,  the  time  at  which  L0  attains  a  shallow  minimum  (see  Figure  5).  Past 
this  time  diffusive  dominance  is  evidenced  by  the  drying  up  of  inner  contours. 


Fig.  3  —  Contour  plots  of  I(x,y)  at  selected  times  illustrating  the  demarcation  between  bifurcating 
(left  column)  and  non-bifurcating  (right  column)  states  for  M  =  10.  Contours  are  spaced  linearly 
between  1  =  1  (ambient)  and  M.  Limits  on  the  critical  R  are  evaluated  at  times  when  bifurcation  of 
the  low  K  striation  (second  plot  in  left  column  for  Figures  1—4)  is  eminent.  Times  appear  in  the 
upper  right  of  each  plot. 
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Fig.  4  —  Contour  plots  of  2{x,y)  at  selected  times  illustrating  the  demarcation  between  bifurcating 
(left  column)  and  non-bifurcating  (right  column)  states  for  M  =  30.  Contours  are  spaced  linearly 
between  Z*1  (ambient)  and  M.  Limits  on  the  critical  R  are  evaluated  at  times  when  bifurcation  of 
the  low  K  striation  (second  plot  in  left  column  for  Figures  1—4)  is  eminent.  Times  appear  in  the 
upper  right  of  each  plot. 


Figure  2.  Curve  1  is  the  bifurcating  case,  and  curve  2  is  the 


Curves  1  and  2  of  Figure  5  shows  L  as  a  function  of  time  for  the  M  -  5 

o 

simulations.  Both  curves  show  an  initial  transient  for  0  £  t  ^  125  sec, 
indicating  that  the  initial  condition  is  far  from  a  "natural"  quasistationary 
state.  For  125  £  t  300  sec,  curve  1  shows  Lq  slowly  decreasing,  correspond¬ 
ing  to  the  well  known  backside  steepening  phenomenon.  Then  for  t  >  300  sec, 

Lq  decreases  rapidly,  indicating  the  onset  of  bifurcation.  This  abrupt  change 

in  the  slope  of  L  (t)  coincides  with  bifurcation  in  all  our  runs.  In  con- 
o 

trast,  curve  2  of  Figure  5  shows  Lq  settling  down  to  a  nearly  constant  value 
after  the  transient  period.  Lq  reaches  a  shallow  minimum  at  300  sec.,  in¬ 
dicating  the  onset  of  diffusive  dominance  rather  than  bifurcation.  The  coin¬ 
cidence  of  this  minimum  with  the  knee  in  curve  1  is  fortuitous,  since  the 
minimum  can  be  made  to  occur  earlier  by  increasing  the  diffusivity. 

The  bounds  on  critical  R  values  gleaned  from  the  simulations  are 
given  in  Table  1.  The  column  denoted  "resolution"  refers  to  results  from 
Figures  1-4  and  a  duplicate  set  of  simulations  performed  with  twice  the 
spatial  resolution  on  the  162  x  82  grid  (boundaries  are  placed  closer  to 
the  structure).  The  degree  to  which  upper  and  lower  bounds  are  different 
for  different  spatial  resolution  gives  a  measure  of  discretization  and 
boundary  effects.  The  geometric  means  of  all  R  values  from  Table  1  are  1207 
and  1355  for  high  and  low  resolution,  respectively,  and  differ  by  only  11%. 
The  geometric  mean  of  upper  and  lower  bounds  for  a  given  M  is  plotted  as  a 
function  of  M  in  Figure  6.  This  plot  suggests  that  the  critical  R  attains 
a  minimum  value  between  600  and  700  for  M  ~  4. 
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Table  1  —  Critical  R  vs.  Conductivity  Ratio 


Resolution 


Range 


Geometric 

Mean 


High 

Low 


1629  > 

R  .  >  993 

2402 

1527 

769 

450 

971 

623 

1110 

692 

1450 

853 

4440 

2310 

2722 

1519 

1561 


676 


987 


2552 


Fig.  6  Critical  R  vs  M.  Solid  points:  mean  values  from  Table  1.  Error  bars  denote  log  variances. 
Circled  points:  density  falloff  scale  size  L  vs  M  from  rocket  data  for  event  Esther  [Table  5.1  of 
Baker  et  al.,  1978] . 


The  shape  of  this  curve  is  a  gratifying  result  which  lends  support  to 
the  simulations.  Linson  [1975J  offered  an  ad  hoc  model  to  explain  observed 
onset  times  as  a  function  of  M.  This  model,  based  on  slip  velocities  for 
elliptical  piecewise  constant  density  ("waterbag")  clouds,  predicted  a  U- 
shaped  curve  with  a  minimum  onset  time  for  some  M  >  2,  the  exact  value  de¬ 
pending  on  the  shape  of  the  cloud.  This  model  was  found  qualitatively  con¬ 
sistent  with  simulation  results  for  realistic  cloud  profiles  [McDonald  et  al. , 
1980],  This  minimum  onset  time  may  be  viewed  as  a  maximum  bifurcation  ten¬ 
dency.  The  amount  of  diffusion  required  to  halt  bifurcation  is  thus  maximum 
at  some  M  value,  resulting  in  a  minimum  critical  R.  A  recent  result  of 
Overman  and  Zabusky  [1980J  for  circular  waterbag  clouds  also  supports  the 
qualitative  dependence  of  the  critical  R  upon  M  as  shown  in  Figure  6.  They 

find  that  shielding  and  dissipation  cooperate  so  as  to  produce  an  effective 

2 

diffusivity  which  is  the  actual  diffusivitv  times  (M+l)  /(M-l).  If  we  assume 

that  the  amount  of  effective  diffusivity  required  to  halt  bifurcation  is  in- 

2 

sensitive  to  M,  the  critical  R  should  be  proportional  to  (M+l)  /(M-l).  This 
expression  has  a  minimum  at  M  =  3,  in  approximate  agreement  with  the  curve 
of  Figure  6.  In  fact,  with  the  exception  of  the  M  =  2  case,  the  expression 
R  =  75  (M+l) /(M-l)  gives  the  mean  R  values  of  Table  1  to  3%  accuracy.  For 
high  M,  the  polarization  charges  responsible  for  shielding  reside  in  a  thin 
layer  in  which  the  cloud's  conductivity  rises  from  ambient  to  a  few  times 
the  ambient  value.  The  thinness  of  this  layer  results  in  agreement  of  the 
cloud's  polarization  electric  field  with  that  of  a  suitable  waterbag  cloud. 
Other  data  from  the  simulations  are  presented  in  Table  2.  The  margin¬ 
ally  stable  structures  in  the  right  columns  of  Figs.  1-4  have  been  Fourier 
analyzed  at  times  at  which  the  critical  R  is  evaluated.  Summing  the  two 
dimensional  spectral  power  over  the  transverse  wavenumber  yields  one  diraen- 
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sional  in  situ  power  spectra  P(k  )  «  k  x  and  P(k  )  a  k  ^  for  the  x  and 

-  xx  y  y 

y  directions  respectively.  The  proportionality  is  valid  only  in  the  central 
portion  of  the  spectrum.  Diffusive  dominance  is  evident  in  the  steepness 
of  the  y  spectra  as  compared  to  the  x  spectra.  This  suggests  the  possibility 
that  a  rocket  moving  transverse  to  the  neutral  wind  through  a  "frozen"  cloud 
might  find  a  power  spectral  density  steeper  by  one  or  two  powers  of  k  than 
a  rocket  moving  parallel  to  the  neutral  wind  (assuming  that  turbulent  fluctu¬ 
ations  can  be  removed  from  the  data).  Finally,  the  coupling  coefficient  Q  is 

given  which  is  simply  the  ratio  of  striation  centroid  velocity  V  to  ambient 

c 

plasma  drift  speed  Vq  in  the  frame  of  the  neutral  atmosphere.  In  agreement 
with  analytic  results  from  a  waterbag  model  ^Linson,  1975],  £  is  roughly 
proportional  to  M  ^  as  a  result  of  electric  field  shielding  at  high  M. 

Figures  1-4  reveal  that  the  manner  in  which  an  unstable  structure 
comes  apart  is  dependent  on  the  conductivity  ratio.  For  M  <  10,  depletions 
are  able  to  penetrate  the  initial  structure  with  ease.  However,  for  M  >  10, 
there  is  a  tendency  for  the  secondary  structure  to  be  confined  to  the  surface 
of  the  original  structure.  High  conductivity  clouds  are  peeled  like  an 
onion,  while  low  conductivity  clouds  are  cut  like  an  apple.  This  effect  is 
even  more  apparent  in  our  earlier  work  [McDonald  et  al.,  1980]. 


Table  2  —  Spectral  Indices  and  Centroid  Drifts  for  Marginally  Stable  States 


M 

TIME 

n 

X 

n 

y 

o 

> 

o 

> 

o 

> 

'•w*' 

^0 

2 

360 

2.2 

5.4 

0.634 

5 

280 

1.8 

4.2 

0.305 

10 

280 

1.8 

3.1 

0.160 

30 

480 

2.0 

3.1 

0.057 
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5.  Summary  and  Comparison  with  Data 


For  plasma  structures  resembling  striation  tips,  the  model  (1)  -  (3) 
implies  that  the  demarcation  between  bifurcating  and  non  bifurcating  states 
is  determined  by  the  value  of  R  (analogous  to  the  Reynolds  number  in  hydro¬ 
dynamics)  and  the  conductivity  ratio  M.  This  demarcation  has  been  estimated 
by  computer  simulation  and  is  presented  in  Figure  6.  The  critical  R  reaches 
a  minimum  of  roughly  700  for  M  ~  4.  Its  dependence  upon  M  is  similar  to 
that  of  the  onset  time  for  striation  emergence  from  a  two  dimensional  cloud 
[Linson,  1975;  McDonald  et  al.,  1980],  This  similarity  is  a  result  of  bi¬ 
furcation  tendency  being  offset  by  a  proportional  amount  of  diffusion. 

Our  results  compare  favorably  with  three  distinct  types  of  experi¬ 
mental  observations.  These  are  a)  minimum  scale  sizes  as  determined  by 
rocket  probe  data;  b)  the  "freezing-up"  of  visible  structure  on  scales  of 
approximately  one  kilometer;  and  c)  the  lifetime  of  the  "frozen"  state. 

We  shall  now  discuss  each  of  these  areas. 


A.  Minimum  Scale  Sizes 

If  M  is  known  for  a  structure,  we  can  find  the  critical  R  from  the 
curve  of  Figure  6.  We  can  then  use  (10)  to  estimate  L^,  providing  we  have 
estimates  for  Vq  and  K.  Since  Lq  is  proportional  to  K,  one  expects  the 
smallest  scales  to  be  determined  by  the  classical  (non  turbulent)  diffusivity 


0  Ve  ckT 
Q  eB  ’ 


(16) 


e 

where  is  the  sum  of  electron  collision  frequencies  with  cloud  ions  and 
ambient  neutrals,  is  the  electron  gyrofrequency,  k  is  Boltzmann's  con¬ 
stant,  T  is  the  plasma  temperature,  and  e  is  the  electron  charge.  The 
expression  (]6)  is  equivalent  to  that  used  in  eq.  (22)  of  Perkins  et  al. 


[1973]. 
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For  plasma  number  densities  n  >  3  x  10  cm  at  altitudes  of  150  km 
or  greater,  is  dominated  by  ionic  collisions.  Thus  K  varies  approximately 

-h  -2 

as  nT  B  ,  as  a  result  of  the  following. 

ve  ~  Vei  =  (34  +  4.18  log10(T3/n))  nT'3/2,  (1?) 

“  3  o  1 

where  n  is  in  cm  ,  T  is  in  °K,  and  is  in  sec  •  High  resolution  plasma 
probe  data  for  the  barium  cloud  Esther  [Figure  2  of  Baker  and  Ulwick,  1978] 

£  CO 

at  approximately  170  km  give  5  x  10  >  n  >  3  x  10  cm  .  We  have  from  (16) 

and  (17)  with  T  =  1000  °K, 


K  =  0.182  m2/s,  n  =  3  x  105  cm-3 

=  2.71  m2/s,  n  =  5  x  106  cm-3  (18) 

In  order  to  estimate  M  for  individual  striations  in  the  cloud,  we  need  £ 

values  in  the  ambient  ionosphere  and  throughout  the  cloud  (denoted  by  £ 

a 

and  £  ,  respectively).  Taking  a  magnetic  dip  angle  I  =  61.5°  for  Eglin 
AFB,  Florida,  Francis  and  Perkins  [1975]  give  for  the  twilight  ionosphere 


£  «  3.5  mho  (19) 

a 

Since  the  cloud  is  contained  well  above  125  km  (where  V.  Hi.  ~  1  for  B  +) 

i  i  a 

we  can  take 


Ci/'cVV1'”1'  <20) 

where  the  path  of  integration  is  along  the  magnetic  field,  h  is  altitude, 
and  field  line  curvature  is  neglected.  Use  of  altitude  rather  than  distance 
along  the  field  facilitates  comparison  with  data.  Assuming 
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2  2 

nc  =  n0  exP  <~h  /D  ). 

=  VQ  exp  (-h/H), 

(20  )  gives 

Ec=D  no  CT  cscI  n  **  exp(D2/4H2).  (21) 

Zero  subscripts  refer  to  values  at  the  cloud  peak  along  a  given  field  line. 

For  conditions  applicable  to  170  km  altitude  at  sunset,  the  result  of 

Linson  and  Baxter  [1977]  for  B  +  gives  V  Hi.  =  .0749  and  H  =  25.0  km. 
-  a  oi 

Taking  B  =  ,5G,  and  D  =  17.31  km  corresponding  to  Esther's  full  width 
at  half  maximum  [Baker  and  Ulwick ,  1978] .of  24  km,  (21)  gives 

E  =  n  (cm  ■*)  x  9.433  x  10  ^  mho 
c  o 

The  measured  values  5  x  10^  i  n  i  3  x  10  cm  give 

47.2  z  Z  2  2.83  mho.  (23) 

c 

M  values  for  the  individual  striations  in  the  Esther  cloud  are  taken  to  be 
ratios  of  E  =  Z  +  E  for  neighboring  peaks  and  valleys  as  given  by  Table 

3  C 

5.1  of  Baker  et  al.  [1978].  Scale  sizes  L  for  density  falloff  are  given  in 
this  table  and  are  plotted  as  circled  points  in  Figure  6.  One  sees  that 

minimum  scale  sizes  for  a  given  M  are  in  agreement  with  our  simulation 
results,  providing  we  take 

K/V  =  0.02  meter.  (24) 

o 

2 

With  Vq  =  100  m/s,  we  have  K  =  2  m  /s,  in  agreement  with  the  range  given 
by  (18).  The  existence  of  scale  sizes  well  above  the  curve  of  Fig.  6  could 
be  a  result  of  oblique  cuts  through  shanks  rather  than  tips  of  striations. 

It  could  also  be  a  result  of  an  enhanced  diffusivity  (to  be  discussed  below),, 
lack  of  complete  temporal  development,  or  effects  not  included  in  the  model 
(as  mentioned  in  Sec.l). 
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B.  Kilometer  Scale  Visible  Structure 
Departing  the  Issue  of  minimum  scale  sizes  which  seem  to  be  controlled 
by  classical  diffusion,  our  results  can  be  used  to  address  the  frequent 
emergence  of  visible  structure  with  scale  sizes  of  order  1  kilometer.  We 
hypothesize  that  plasma  turbulence  may  be  involved  in  these  structures.  This 
admits  the  possibility  of  individual  parcels  of  plasma  macroscopically  acting 
as  agents  of  fluid  transport,  while  microscopically  being  subject  to  classical 
particle  diffusion.  Thus  the  smallest  structures  may  reflect  electron 
diffusion  dominance,  while  much  larger  structures  may  reflect  enhanced 
diffusivity  effects.  In  order  to  model  these  large  scale  structures,  we 
replace  the  classical  particle  diffusivity  (16)  by  an  enhanced  diffusion  term 
descriptive  of  a  turbulent  process.  Our  goal  is  not  to  prove  that  turbulence 
is  operative,  but  to  offer  at  least  one  physically  plausible  explanation  for 
observed  scale  sizes  and  lifetimes. 


If  we  take  Bohm  diffusion  [Chen,  1974] 


i_  ckT 
16  eB 


(25) 


as  descriptive  of  turbulent  diffusivity,  we  find  K  values  much  greater  than 

those  given  by  (16)  for  typical  barium  releases.  With  T  =  1000°K  and 

2 

B  »  .  5  G,  (25)  gives  Kb=  108  m  / s.  The  curve  in  Figure  6  suggests  that 

R  «  1000  is  accurate  to  +  30%  for  a  range  of  M  values  from  approximately  2  to 

15.  Thus  (9)  would  predict  Lq  «  1  km  for  the  gradient  scale  size  for  a 
variety  of  typical  barium  releases  in  which  e*  100  m/s.  This  in  agreement 

with  observations  cited  in  the  Introduction  [J.  A.  Fedder,  W.  Chestnut, 

private  communication,  1980]. 


26 


C.  Lifetime  of  the  Kilometer  Scale  Structure 


If  the  persistent  structure  is  viewed  as  marginally  stable  to  bifur¬ 
cation,  its  lifetime  will  be  determined  by  diffusion,  and  should  be  in 

2 

approximate  agreement  with  (11).  Using  Lq  =  1km  and  kb  =  100  to  /s,  we  find 

4 

a  diffusion  time  scale  of  10  sec.  =  2.7  hours.  Horizontal  drift  and  sun¬ 
lighting  conditions  make  visual  tracking  of. a  cloud  for  this  length  of  time 
very  difficult.  However,  satellite  communications  experiments  have  been 

conducted  for  timescales  on  this  order  during  the  STRESS  program.  For  the 
five  STRESS  releases,  propagation  through  the  clouds  resulted  in  signal 

fading  for  durations  of  1.4  to  3.7  hours  after  release  [Prettie  et  al.  1977]. 

The  beam  frequency  of  341  MHz  and  the  approximate  200  km  range  from  cloud  to 

receiving  aircraft  implies  that  the  fading  was  sensitive  to  structure  on 

u 

scales  roughly  that  of  the  Fresnel  size  (Xz)^  ■  420  meters. 
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OICY  ATrN  XPFS  MAJ  3.  SrEPHAN 
OICY  ATtn  ADWATE  MAJ  BRUCE  BAUER 
OICY  ArTN  NRT 

OICY  ArTN  OOK  CHIEF  SC  I £** I S T 

SAMSO/SK 
P.  O.  BOX  92960 
rfORLXAY  POS’Au  CENTER 
.OS  ANGElES,  CA  90009 

OICY  A*rN  SKA  'SPACE  COMM  SYSTEMS)  M.  CuAVlN 
SAMSO/MN 

NORTON  AFB,  CA  92*09 
^(NUrEMAN) 

OICY  A*TN  LrC  KENEDY 


COMMAMOER 

ROME  AJR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  EEP  A.  LORENTZEN 

department  of  energy 
ALBUQUERQUE  OPERATIQNS  OFFICE 
p.  O.  BOX  5400 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  D.  SHERWQOO 

DEPARtMENT  OF  ENERGY 
LIBRARY  ROOM  G-042 
WASHINGtON,  D.C.  20545 

OICY  AT^n  DOC  CON  FOR  A.  LA80WITZ 

EG&G,  INC. 

LOS  ALAMOS  DIVISION 

P.  O.  BOX  809 

LOS  ALAMOS,  NM  85544 

OICY  ATTN  DOC  CON  FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.  0.  BOX  808 
LIVERMORE,  CA  94550 

OICY  AtTn  DOC  CON  FOR  TECH  INFO  DEPT 
OICY  ATTN  DOC  CON  FOR  L-389  R.  0TT 
OICY  ATTN  DOC  CON  FOR  L-31  R .  HAGER 
OICY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 


LOS  ALAMOS  SCIENTIFIC  LABORATORY 
P.  0.  BOX  1663 


AlMOS 

,  ^ 

87545 

OICY 

ATTN 

DOC 

CON 

FOR 

J, 

WOLCOTT 

OICY 

ATTN 

DOC 

CON 

FOR 

R. 

F.  TASCHEK 

OICY 

ATTN 

DOC 

CON 

FOR 

E. 

JOFCS 

OICY 

attn 

DOC 

CON 

FOR 

U. 

FVALIK 

OICY 

attn 

DOC 

CON 

FOR 

R. 

JEFFRIES 

OICY 

ATTN 

DOC 

CON 

FOR 

J. 

Z1MW 

OICY 

ATTN 

OOC  CON 

FOR 

e. 

KEATON 

flier 

ATT*  DOC 

CON 

FOR 

D. 

WESTERVELT 

SANDIA  LABORATORIES 
P.  0.  BOX  5800 

Albuquerque,  nm  87115 

OICY  ATTN  DOC  CON  FOR  J.  MARTIN 

OICY  ATTN  DOC  CON  FOR  W.  BROWN 

OICY  ATtn  DOC  CON  FOR  A,  TH0RN8R0UGH 

OICY  ATtN  DOC  CON  FOR  T.  WRIGHT 

OICY  AT’’*  DOC  CON  FOR  D.  DAHLGREN 

OICY  AT^N  DOC  CON  FOR  3141 

OICY  Attn  ooc  CON  FOR  SPACE  PR0JECT  DlV 

SAND I A  LABORATORIES 
LIVERMORE  uABORATORY 
P.  0.  BOX  969 
LIVERMORE,  CA  94550 

OICY  A^N  DOC  CON  FOR  B.  MURPHEY 
OICY  A"N  DOC  CON  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
department  OF  ENERGY 
WASHINGTON,  D.C.  20545 

OICY  ATtn  DOC  CON  FOR  D.  GALE 


OTHER  GOVERNMENT 


CENTRAL  INTELLIGENCE  AGENCY 
A’”’*  RD/SI,  RM  5G48,  mQ  BLOG 
WASHINGTON,  D.C.  20505 

OICY  ATTN  OS  I /PS  ID  RM  5F  19 

OEPART^NT  OF  COMf€RCE 
NATlONAL  BUREAU  OF  STAJ*>ARDS 
wAShINGtON,  D.C.  20234 

(ALL  CORRES:  at^n  SEC  Officer  for) 
OICY  a**'*  R.  MOORE 
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IN$T1TUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  S  InFQ  AOMIN 
BOULDER,  CO  80503 

OICY  ATTN  A.  JEAN  (UNCLASS  Of*_Y) 
OICY  ATTN  W.  U'T-AUT 
01CY  ATTN  D.  CROMBIE 
01CT  ATTN  L.  BERRY 

NAT IONAL  OCEANIC  £  ATMOSPHERIC  AOMIN 
ENVIRO<^€NTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 

01CY  ATTN  R.  GRUBB 

01CY  ATTN  AERONOMY  LAS  G.  REID 


DEPARTMENT  of  defense  contractors 


AEROSPACE  corporation 


0.  BOX  92957 
ANGELES,  CA 

90009 

OICY 

ATTN 

I.  GARFUMCEL 

OICY 

ATTN 

T.  salmi 

OICY 

ATTN 

V.  JOSEPHSON 

OICY 

ATTN 

S.  BOWER 

OICY 

ATTN 

N.  STOCKWELL 

OICY 

ATTN 

0.  OLSEN 

OICY 

ATTN 

J.  CARTER 

OICY 

at^n 

F.  MORSE 

OICY 

AT^N 

SMFA  FOR  FWW 

ANALYTICAL  systems  engineering  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

01CY  ftADio  SCIENCES 


BERKELEY  RESEARCH  ASSOCIATES,  INC. 

P.  0.  BOX  983 
BERKELEY,  CA  94701 

01CY  ATtn  J.  WORKMAN 

BOEING  COMPANY,  THE 
P.  0.  BOX  3707 
SEAttLe ,  wA  98124 

OiCY  At^  G.  KEISTER 
01CY  ATTN  0.  MURRAY 
OICY  AtTn  G.  HALL 
OICY  ATTN  J.  KENNEY 

CALIFORNIA  AT  SAN  DIEGO,  UNIV  OF 
I PAPS,  0-019 
LA  JOLLA,  CA  92093 

OICY  ATtn  HENRY  G.  BOOKER 

BROWN  ENGINEERING  COMPANY,  INC. 
CURINGS  RESEARCH  PARK 
HUNTSVILLE,  AL  35807 

OICY  ATTN  ROMEO  A.  DELISERIS 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  technology  SQUARE 
CAMBRIDGE,  MA  02139 

OICY  ATTN  D.  3.  COX 
OICY  ATTN  J.  P.  GILMORE 

computer  SCIENCES  CORPORATION 
6565  ARL[NGTON  0LVD 
FALLS  CHURCH,  VA  22046 
OICY  ATTN  M.  BLANC 
OICY  ATTN  JOHN  SPOOR 
OICY  ATTN  C.  NAIL 

COMSAT  LABORATORIES 

linthicum  rqao 
CLARKSBURG,  MD  20734 
OICY  ATTN  G.  HYDE 

CORNELL  UNI  VERS  I TY 

department  of  electrical  engineering 

ITHACA,  NY  14850 

OICY  ATTN  D.  T.  FARLEY  JR 


electrospace  systems,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 

OICY  AT’rN  M.  _OGSTON 

OICY  A^N  SECURITY  (PAUL  philips; 

ESL  INC. 

495  JAVA  DRIVE 
Sunnyvale,  ca  94086 

OICY  ATTN  J.  RO0ERTS 
OICY  A^N  JAMES  MARSHALL. 

OICY  AT^N  C.  W.  PR£T'r l£ 

FORD  AEROSPACE  S  COMNJNICA’ 1 ONS  CORP 
3939  FABIAN  -AY 
PALO  ALTO,  CA  94303 

OICY  AT^n  j.  T.  MA^IN&lEY 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GOOOARD  BLVD  KING  OF  PRUSSIA 
P.  0.  BOX  3555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.  h.  BORtNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.  0.  BOX  1122 
SYRACUSE,  NY  13201 

OICY  A^N  F.  RE  I  BERT 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  SADIES 
816  state  stree^  (P.O.  DRAWER  QQ) 

SANTA  BARBARA,  CA  95102 
OICY  ATtN  DAS  I  AC 
OICY  ATTN  DON  CHANDLER 
OICY  A^n  tom  BARRE^ 

OICY  ATTN  TIM  STEPHANS 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATtn  8.  GAMBI-L 
OICY  ATTN  MACK  SrANT0N 

GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC. 

HM6S 

COURT  STREET 
SYRACUSE,  NY  13201 

OiCY  ATTN  G.  MILLMAN 

GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
P.  0.  BOX  6770 
SANTA  BARBARA,  CA  93111 

OICY  ATTN  JOHN  ISE  JR 
OICY  ATTN  JOEL  GARBARINO 


GEOPHYSICAL  iNSTinjTE 
UNIVERSrY  OF  ALASKA 
FAIRBAMCS,  AK  99701 

(ALL  CLASS  AtTN:  SECURITY  OFFICER) 
OICY  ATTN  T.  N.  DAVIS  (UNCL  ONLY) 

OICY  ATTN  NEAL  BROW  (UNCL  ONLY) 

OICY  ATTN  TECHNICAL  LIBRARY 

GTE  SYL VANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP -EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

OICY  ATTN  MARSHAL  CROSS 

ILLINOIS,  UNIVERSITY  OF 
DEPARThCNT  OF  ELECTRICAL  ENGINEERING 
URSANA,  (L  61803 

OICY  ATTN  K.  YEH 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
801  S.  WRIGHT  STREET 
UR8ANA,  IL  60680 

(All  CORRES  ATTN  SECURITY  supervisor  for) 
OICY  ATTN  K.  YEH 
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INS’l'U’E  FOR  DEFENSE  ANALYSES 
400  ARMY -NAVY  DRIVE 
ARLINGTON,  VA  22202 

01CY  ATTN  J.  M.  AE1N 
01CY  ATTN  ERNEST  SAUER 
OKY  AT’N  MANS  WOLFHARD 
OICY  ATTN  JOEL  8ENGSTON 

HSS,  1NL. 

2  Alfred  circle 

BEDFORD,  MA  01750 

OICY  A’TN  DONALD  MANSEN 


INTL  TEL  S  TELEGRAPH  CORPORATION 

500  washing’on  avenue 

nu’lEy,  nj  07110 

OICY  AT’-n  technical  library 


JAYCOR 

1401  CAMINO  del  MAR 
DEL  MAR,  CA  92014 

OICY  at’n  S.  r.  Goldman 


JOHNS  HOPKINS  'JN I  VERS  I  TY 

applied  physics  laboratory 

JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
OICY  ATTN  THOMAS  potemra 
OICY  Attn  JOHN  DASSOULAS 

lOCKHEEO  missiles  S  SPACE  CO  INC 
p.  O.  BOX  504 
SUNNYVALE,  CA  94088 

OICY  AT’N  DEPT  60-12 
OICY  ATTN  D.  R.  CHURCHILL 

LOCKHEEO  MISSILES  AND  SPACE  CO  INC 
5251  HANOVER  STREET 
PALO  ALTO,  CA  94504 

OICY  AT’n  MARTIN  WALT  DEPT  52-10 
OICY  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
OICY  ATTn  W.  L.  IIHOF  DEPT  52-12 

KAMAN  SCIENCES  CORP 
P.  0.  BOX  7465 
COLORAOO  SPRINGS,  CO  80955 
OICY  A”N  T.  MEAGHER 

LIF*CABIT  CORP 
10455  ROSELLE 
SAN  DIEGO,  CA  92121 

OICY  A”N  IRWIN  JACOBS 

M.I.T.  LINCOLN  LABORATORY 
P.  0.  BOX  75 
LEXINGTON,  MA  02175 

OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  P.  WALDRON 
OICY  ATTN  L.  LOUGHLIN 
OICY  ATTN  D.  CLARK 

MARTIN  MARIET’A  CORP 
ORLANDO  DIVISION 
P.  0.  BOX  5857 
ORLANOO,  fl  52805 

OICY  ATTN  R.  HEFFNER 

MCOOM4ELL  OOUGLAS  CORPORATION 
5501  BOLSA  AVENUE 
HWTINGTON  BEACH,  CA  92647 
OICY  ATTN  N.  HARRIS 
OICY  ATTN  J.  MOULE 
OICY  ATTN  GEORGE  MROZ 
OICY  AT’N  W.  OLSON 
OICY  AT’N  R.  W.  HALPRIN 
OICY  ATTN  TECHNICAL  LIBRARY  SERVICES 


MISSION  RESEARCH  CORPORATION 


735  STATE 

S’REET 

SAN’ A  BARBARA,  CA  93101 

OICY 

ATTn  P.  FISCHER 

OICY 

ATTN  W.  F.  CREVIER 

OICY 

ATTN  STEVEN  L.  GUTSCHE 

OICY 

attn  d.  sappenf i eld 

OICY 

attn  R.  BOGUSCH 

OICY 

attn  R.  HENDRICK 

OICY 

attn  RALPH  <ILB 

OICY 

AT’n  DAVE  SOWLE 

OICY 

ATTN  f.  FAJEN 

OICY 

ATTn  m.  SCHE18E 

OICY 

ATTN  CONRAD  L.  LONGM1RE 

OICY 

ATTn  WARREN  a.  schlueter 

M!TR e  CORPORATION,  THE 

P.  0.  BOX 

208 

BEDFORD,  * A  01730 

OICY 

attn  JOHN  .MORGANS  TERN 

OICY 

ATTN  G.  HARDING 

OICY 

A’Tn  C.  E.  CALLAHAN 

Ml ’RE  CORP 

WES’" GATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 

OICY  ATTN  W.  HALL 
OICY  ATTN  W.  FOS’ER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERFIELD  BLVD. 

SANTA  MONICA,  CA  90404 

OICY  AT’N  £.  C.  FIELD  JR 

PENNSYLVANIA  S’ATE  UNlVERSITY 
IONOSPHERE  RESEARCH  LAB 
518  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 

(NO  CLASSIFIED  TO  ’HIS  ADDRESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PhOTomE’RICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02175 

OICY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS  INC. 
p.  0.  BOX  5027 
BELLEVUE,  WA  98009 

OICY  A’Tn  E.  J.  FREMOUW 

PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  10367 
OAKLAND,  CA.  94610 

attn:  a.  THOMSON 


R  S  D  ASSOCIATES 

P.  0.  BOX  9695 

MARINA  DEL  REY,  CA  90291 

OICY  ATTN  FORRES’  GILMORE 
OICY  ATTN  BRYAN  GABBARD 
OICY  ATTN  WILLIAM  B.  WRIGHT  JR 
OICY  ATTN  ROBERT  F.  LELEVIER 
OICY  ATTN  WILLIAM  J.  KAR2AS 
OICY  AT’N  H.  ORY 
OICY  ATTN  C.  MACDONALD 
OICY  ATTN  R.  TURCO 

RANO  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  8EDR02IAN 

RIVERSIDE  RESEARCH  INSTI’uTE 
80  WEST  ETC  AVENUE 
NEW  YORK,  NY  10025 

OICY  AT’N  VINCE  TRAPANI 
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SCIENCE  APPLICATIONS,  INC. 

P.  0.  BOX  2351 
lA  JOLLA,  CA  92038 

OICY  ATTN  LEWIS  M.  LINSON 
01CY  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  0.  SACHS 
01CY  ATTN  E.  A.  STRAKER 
01CY  ATTN  CURTIS  A.  SMITH 
01CY  ATTN  JACK  MCDOUGALL 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01/76 

01CY  ATTN  BARBARA  ADAMS 

Science  Applications,  Inc 
1700  Goodridge  Dr. 

McLean,  Va.  22102 

01  cy  ATTN  J.  COCKAYNE 


SRI  INTERNATIONAL 
333  RAVENSWOOO  AVENUE 
MENLO  PARK,  CA  94025 

01CY  ATtn  DONALD  NElLSON 

oicy  attn  alan  8urns 
01CY  ATTN  g.  smith 
OICY  ATTN  l.  L.  COBB 
OICY  attn  DAVID  A.  JOHNSON 
OICY  ATTN  WALtER  G.  CHESNU’’ 
OICY  ATTN  CHARLES  L.  RINO 
OICY  ATTN  WALTER  JAYE 
OICY  ATTN  M.  BARON 
OICY  attn  RAY  L.  leaoabrano 
OICY  ATTN  G.  CARPENTER 
OICY  AT^N  G.  PRICE 
OICY  ATTN  J.  PETERSON 
OICY  ATTN  R.  HAKE,  JR. 

OICY  ATTN  V.  GONZALES 
OICY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OICY  AT^n  W.  P.  80QUI St 


Lockheed  Missile  &  Space  Co.,  Inc. 
Huntsville  Research  &  Engineering 
Center 

4800  Bradford  Dr. 

Huntsville,  Ala.  35807 
01  cy  ATTN  DALE  H.  DIVIS 


TRW  DEFENSE  £  SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONDO  BEACH,  CA  90278 

QICY  ATTN  R .  < .  PLEBUCH 
OICY  AT-n  S.  ALTSCHULER 
OICY  AT^n  0.  DEE 

VISIDYNE,  INC. 

19  THIRD  AVENUE 

NORtM  WES'r  INDUSTRIAL  PARK 

BURLINGTON,  MA  01803 

oicy  at-^n  Charles  Humphrey 
oicy  a-^n  j.  w.  carpenter 
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UNCLASSIFIED  ONLY 


PLEASE  DISTRIBUTE  ONE  COPY  TQ  EACH  OF  ThE  FOLLOWING  PEOPLE: 
ADVANCED  RESEARCH  PROJECTS  AGENCY  (ARPA) 

strategic  technology  office 

ARLINGTON,  VIRGINIA 

CAPT.  DONALD  M.  LEVINE 


COMMAPOER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  of  The  NAVY 
WASHINGTON,  D.C.  20560 

DR.  T.  CZUBA 


NAVAL  RESEARCH  LABORATORY 

WASHINGTON,  D.C.  20375  HARVARD  UNIVERSITY 

HARVARD  SQUARE 

DR.  P.  MANGE  CAMBRIDGE,  MASS.  02138 

DR.  R.  MEIER 

DR.  E.  S2USZCZEWICZ  -  COOE  4127  D*-  M*  B-  MCELROY 

DR.  R.  LIPOZEN 


DR.  o.  GOODMAN  -  COOE  7560 


SCIENCE  APPLICATIONS,  INC. 
1250  prospect  plaza 
la  JOLLA,  CALIFORNIA  92037 

DR.  D.  A.  mAMLIN 
OR.  L.  u { NSON 
DR.  D.  SACHS 


DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL  LABORATORY 
NOAA 

BOULDER,  COLORADO  80302 

DR.  A.  GLENN  JEAN 
DR.  G.  W.  ADAMS 
DR.  D.  N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.  F.  DONNELLY 


PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PENNSYLVANIA  16802 


DR.  J. 

S.  NISBET 

DR. 

P. 

R.  ROHRBAUGH 

DR. 

D. 

E.  8ARAN 

DR. 

L. 

A.  CARPENTER 

DR. 

M. 

LEE 

DR. 

R. 

DIVANY 

DR.  P. 

BENNETT 

DR. 

E. 

KL EVANS 

UNIVERSITY  OF  CALIFORNIA,  LOS  ANGELES 

405  HtLLGARD  AVENUE 

LOS  ANGELES,  CALIFORNIA  90024 

DR.  F.  V.  CORONITI 
DR.  C.  KENNEL 


A.  f.  GEOPHYSICS  LABORATORY 
l.  G.  HANSOM  FIELD 
BEDFORD,  MASS.  01730 


UNIVERSITY  OF  CALIFORNIA,  BERKELEY 
BERKELEY,  CALIFORNIA  94720 

DR.  M.  HUDSON 


DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  «.  SAGALYN 
DR.  J.  M.  FORBES 
DR.  T.  j.  KENESHEA 
DR.  u.  AARONS 


OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VIRGINIA  22217 

DR.  H.  MULLANEY 


UTAH  STATE  UNIVERSITY 
4TH  N.  AND  8TH  STREETS 
LOGAN,  UTAH  84322 

DR.  P.  M.  BANKS 
DR.  R.  HARRIS 
DR.  V.  PETERSON 
DR.  R.  MEGILL 
DR.  K.  BAKER 


CORNELL  UNIVERSITY 
ITHACA,  NEW  YORK  14850 


COMMANDER 

NAVAL  ELECTRONICS  LABORATORY  CENTER 
SAN  DIEGO,  CALIFORNIA  92152 

DR.  M.  BLEIWEISS 
DR.  I.  ROTHMULLER 
DR.  V.  HILDEBRAND 
MR.  R.  ROSE 


DR.  W.  E.  SWARTZ 
DR.  R.  SUDAN 
DR.  D.  FARLEY 
DR.  M.  KELLEY 


NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREEN8ELT,  MARYLAPD  20771 


U.  S.  ARMY  ABERDEEN  RESEARCH  AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  .MARYLAND 


DR.  S.  CHANDRA 
DR.  K.  MAEDO 


PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  NEW  JERSEY  08540 

DR.  P.  PERKINS 
DR.  E.  FRIEMAN 


institute  for  defense  analysis 

400  ARMY /NAVY  DRIVE 
ARLINGTON,  VIRGINIA  22202 

DR.  E.  BAUER 


UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20742 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA.  15215 

DR.  N.  ZA8USKY 
OR.  M.  BlOfDI 


UNIVERSITY  OF  CALIFORNIA 

LOS  ALAMOS  SCIENTIFIC  LABORATORY 

J-10,  MS-664 

LOS  ALAMOS,  NEW  MEXICO  87545 

M.  P0NGRAT2 
0.  SIMONS 
G.  BARASCH 
L.  DUNCAN 
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